Direct measurement of the van der Waals interaction between two single atoms 
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We report on the direct measurement of the van der Waals interaction between two isolated, 
single Rydberg atoms separated by a controlled distance of a few micrometers. By working in 
a regime where the single-atom Rabi frequency of the laser used for excitation to the Rydberg 
state is comparable to the interaction energy, we observe a partial Rydberg blockade, whereby the 
time-dependent populations of the various two-atom states exhibit coherent oscillations with several 
frequencies. A quantitative comparison of the data with a simple model based on the optical Bloch 
equations allows us to extract the van der Waals energy, and to observe its characteristic Ce/R^ 
dependence. The magnitude of the measured Ce coefficient agrees well with an ab-mitio theoretical 
calculation, and we observe its dramatic increase with the principal quantum number n of the 
Rydberg state. Our results not only allow to test an important physical law, but also demonstrate 
a degree of experimental control which opens new perspectives in quantum information processing 
and quantum simulation using long-range interactions between the atoms. 



The van der Waals-London interaction C/vdw = Ce /R^ 
between two neutral, polarizable particles separated by a 
distance R is ubiquitous in nature. It is at the origin of 
many effects in physics, chemistry and biology, ranging 
from the condensation of non-polar gases into liquids, to 
the adhesion properties of gecko toes [1]. Measurements 
of the van der Waals interaction between an atom and a 
metallic or dielectric surface have been carried out with 
a variety of techniques [2], some of which make use of 
Rydberg atoms (atoms with a large principal quantum 
number n), that exhibit very strong interactions [3, 4]. 
However, in this geometry, the sum over the underlying 
Cg / interactions between individual particles results in 
a potential scaling as l/L'^, where L is the atom-surface 
distance. To date, no direct measurement of the van 
der Waals interaction between two isolated, single atoms 
separated by a controlled distance, has been performed. 
Such a measurement would be extremely challenging for 
ground-state atoms due to their very small interaction 
strength. Atoms excited to Rydberg states show dramat- 
ically enhanced interactions, yielding sizable effects even 
for atom separations as large as several micrometers. 

Here, we report on the measurement of the Cg /R^ in- 
teraction between two isolated Rydberg atoms prepared 
in a well defined quantum state. For that purpose, we ir- 
radiate the pair of atoms with a laser coupling the ground 
state \g) and the targeted Rydberg state \r) with a Rabi 
frequency Q.. Depending on the relative strength of C/vdw 
and Ml, two limiting cases can be identified. In the limit 
UvdW ^ fi^j the two atoms behave independently and 
the doubly excited state \rr) can be populated with unit 
probability. On the contrary, if the laser coupling hfl 
is much smaller than J7vdW) the excitation of \rr) is off- 
resonant and thus suppressed (Fig. la), giving rise to the 
so-called Rydberg blockade mechanism [5-7]. This effect 



leads to the appearance of a "blockade sphere" inside of 
which only a single Rydberg excitation can be created, 
giving a simple picture for non-trivial many-body effects 
in driven systems. The Rydberg blockade has been ob- 
served in recent years in extended ensembles [8-12] as 
well as between two atoms [13, 14]. 

When working in the transition region hn ^ C/vdWj 
i.e. in a regime of partial blockade, the blockade sphere 
picture is too simplistic: the populations of the various 
many-atom states undergo coherent collective oscillations 
with several frequencies which depend on L^vdw- In our 
two-atom experiment, this allows us to extract the inter- 
action strength from their time dependence. Our results 
prove that a detailed experimental control over the inter- 
actions between Rydberg atoms is possible. Such control 
is a fundamental prerequisite in view of future applica- 
tions to high-fidelity quantum information processing [15] 
and quantum simulation [16]. 

We use two single ^''Rb atoms at a temperature of 
about 50 /iK, loaded in 1 mK-deep microscopic optical 
dipole traps from a magneto-optical trap (Fig. lb) [17]. 
Our new experimental setup is designed to overcome the 
limitations of the apparatus used in our early studies 
of Rydberg blockade [14] and entanglement [18]. We 
use an aspheric lens [19] with a numerical aperture of 
0.5 (focal length 10 mm, working distance 7 mm) to fo- 
cus two 850 nm trapping beams down to 1.1 /^m (1/e^ 
radius). The distance R between the two traps can 
be varied by changing the incidence angle of the trap- 
ping beams on the lens [20]. The same aspheric lens 
is used to collect the atom fiuorescence from each trap 
on separate avalanche photodiodes (APDs) operating in 
the single-photon counting mode. This allows us to de- 
tect the presence or absence of an atom in each trap. 
To minimize stray electric fields detrimental to Rydberg 
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FIG. 1: a: Two-atom level scheme underlying the Rydberg blockade mechanism. The single-atom Rabi coupling n between \g) 
and |r) gives rise to a collective coupling \/2Q, between \gg), {\gr) + \rg))/\/2, and \rr). b: Sketch of the experimental setup. 
Two single atoms are trapped in microscopic optical dipole traps separated by a distance i? of a few micrometers. Electric field 
control is provided by a set of eight electrodes (four of which, facing the ones displayed here, are not shown). Inset: details of 
the two-photon excitation scheme from \g) to \r). The intermediate-state detuning is A ~ 27r x 740 MHz. 



state manipulations, the surface of the lens facing the 
atoms is coated with an indium-tin oxide layer. In addi- 
tion, eight independent electrodes, similar to the ones 
described in Ref. [21], allow us to cancel stray elec- 
tric fields (which could be as high as 0.15 V/cm). We 
found this cancellation to be critical to the success of the 
experiment. The two-photon excitation sequence from 
\g) — 1551/2, F = 2,tof = 2) to the Rydberg state 
|r) = \nD2/2,mj = 3/2) (see inset of Fig. lb), described 
in [22], allows us to perform coherent oscillations with a 
single-atom Rabi frequency f2/(27r) in the range 500 kHz- 
5 MHz [23], as shown in Fig. 2 A (green triangles). The 
two-photon excitation pulse has a duration r, and the 
laser frequencies are adjusted so that the (light-shift cor- 
rected) one-atom detuning (5 ~ 0. After the excitation 
sequence, we infer the state of each atom by detecting 
its presence or absence in the trap (an atom in the Ry- 
dberg state is slightly anti-trapped by the optical po- 
tential, which results in the loss of the atom). We thus 
reconstruct the populations Pij of the two-atom states 
\ij) {i and j taking the values g and r), by repeating 
each sequence for about 100 times [22]. 

Figure 2(b-e) shows the probabilities P^g + Pgr to ex- 
cite one and only one atom (left column, blue diamonds) 
and Prr to excite both atoms (right column, black disks) 
as a function of the area fir of the excitation pulse, for 
various choices of R and ^l, in the case n = 62. In row 
b, the atoms are far apart (i? ~ 15 /im) and thus almost 
independent. Indeed, the joint probability Pj-r of excit- 
ing both atoms is nearly equal to the product Pi,r ■ P2.r 
(red squares of row a) , where Pi^r is the probability for 
atom i to be in jr), obtained in a control experiment 
where only trap i is present. In this case Prg + Pgr is 
expected to oscillate sinusoidally at frequency 2n be- 
tween and 1/2, close to what we observe experimen- 
tally. On the contrary, in row e, the atoms are close 



to each other {R ~ 4.0 fim) and the Rabi frequency 
(fi ~ 27r X 1.9 MHz) is small enough for almost per- 
fect blockade to occur: at all times, Prr is negligible 
{Prr < 0.06, with an average of 0.036), and thus differs 
substantially from the product of the single-atom exci- 
tation probabilities expected for independent atoms. At 
the same time, Prg + Pgr oscillates at (1.49±0. 07)51, close 
to the collective Rabi frequency y^fi. 

In between those two extreme regimes (rows c and d), 
the blockade is only partial: both Prg + Pgr and Prr 
show a more complex behavior, revealing a non-trivial 
two-atom state. Due to the interactions, the populations 
oscillate in time with several frequencies. We model the 
evolution of the system by solving the optical Bloch equa- 
tions (OBEs) for the two atoms [20, 24, 25]. The single- 
atom Rabi frequency and the (small) damping coef- 
ficient 7 (essentially due to spontaneous emission from 
the 5Pi/2 intermediate state) appearing in the model are 
obtained by fitting the single-atom Rabi oscillation data 
(green triangles in Fig. 2a) . The only remaining parame- 
ters in the model are the van der Waals interaction shift 
Uvdw and the detuning S. We treat C/vdw as an ad- 
justable parameter to fit the solution of the OBEs to 
the two-atom data Prg{t) + Pgr{t) and Prr{t). Examples 
of such fits are presented as solid lines in Fig.2b-e (the 
shaded area shows the confidence interval in the fitted 
value of f/vdw)i and the results are summarized in the 
table of Fig. 2f. We also treat (5 as a free parameter, in 
order to account for slow drifts of the laser system [26] . 
With this method, we obtain the magnitude [J/vdwl: as 
for (5 ~ the sign of the interaction shift does not af- 
fect the time dependence of the populations. We have 
checked that if we deliberately set i5 = and 7 = in the 
fitting procedure (thus reducing our analysis to simply 
solving an effectively three-level Schrodinger equation), 
the values of [J/vdwl obtained depart by at most 20% 
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FIG. 2: From almost independent to fully blockaded atoms, via the partial blockade regime, for n — 62. a: single-atom Rabi 
oscillation Pi,r(f2r) between \g) and \r) (left panel, green triangles), and product of such excitation probabilities for single 
atoms trapped in traps 1 and 2 while the other trap is off (right panel, red squares), b— e: Probability Prg + Pgr to excite one 
and only one atom to the Rydberg state (left panels, blue diamonds), and double Rydberg excitation probability P^r (right 
panels, black disks), versus the excitation pulse area fir. From b to e, the distance R between the atoms decreases, yielding 
increasing blockade. The solid lines show the result of fits of the data by the solution of the optical Bloch equations discussed 
in the text. The shaded areas correspond to one standard deviation in the determination of If/vdwl- f: parameters and results 
of the fits for the data in b— e. The statistical error bars on the determinations of are at the 1% level. 



from the ones obtained with the finite value of 7 and 
with S variable. We have also checked that the interac- 
tion energy yielded by the fits does not depend on the 
chosen Rabi frequency by doubling or halving for sev- 
eral values of R. We emphasize that the convergence of 
the fit is optimal when Uvdw ~ Combined with our 
range of accessible values for fi, this means that in prac- 
tice we can determine values of |J7vdwl/^ in the range 
0.1-10 MHz. 

Figure 3 is a double-logarithmic plot of the values 
of |?7vdw| extracted from such fits versus the inter- 
atomic distance R, for three different Rydberg states 
|r) = |nD3/2,'7ij = 3/2) with principal quantum num- 
bers n = 53, n = 62 and n = 82. The data is con- 
sistent with a power-law behavior of exponent 6. Our 



results illustrate the dramatic dependence of the interac- 
tion strength with the principal quantum number n: for 
instance, changing n from 53 to 62 at a given R yields a 
50-fold increase in the interaction. 

The data also show good agreement with the theoreti- 
cal expressions of | J7vdw I (solid lines in Fig. 3) that we ob- 
tained by direct numerical diagonalization of the hamilto- 
nian of two atoms interacting via the dipole-dipole inter- 
action (DDI) [20]. At the large distances relevant here, 
the effect of the DDI is simply to shift the two-atom 
Rydberg levels by a quantity Cg/i?^, that can be ob- 
tained from second-order perturbation theory. At lower 
distances, mixing between adjacent Rydberg levels would 
occur [27], altering the 1/R^ character of the interaction 
(this can already be seen on the curve for n = 53 for 
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FIG. 3: Measurement of the van der Waals interaction. The 
magnitude |(7vdw| of the interaction energy between two 
atoms in the \nD2,/2,nn,j = 3/2) state, as extracted from the 
fits of Fig. 2, is plotted as a function of the distance R be- 
tween the atoms, for n = 53 (blue), n = 62 (red), and n = 82 
(green). The vertical error bars are one standard deviation 
in the fitted value of [(yvdwl- The solid lines correspond to 
the ab-imtio theoretical expectation (see text) without any 
adjustable parameter; the shaded areas illustrate the 5% sys- 
tematic uncertainty in our calibration of R. 



i? < 4 /im). We have checked that the values of the Cq 
coefficients that we obtain in our calculation are in good 
agreement with the ones reported in [28]. 

In summary, using coherent oscillations in the four- 
level system constituted by two atoms, we have directly 
measured the van der Waals interaction energy between 
two isolated Rydberg atoms. The level of control demon- 
strated here opens exciting perspectives in multi-atom 
systems, for the observation of geometry-dependent ef- 
fects due to the anisotropy of the dipolar interaction [29] 
or the non-additivity of van der Waals interactions [30] . 
It is also a prerequisite to the application of optimal 
strategies for generating high-fidelity, many-atom entan- 
glement using the Rydberg blockade, as well as for quan- 
tum simulation of long-range interacting spin systems. 
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Supplementary material 



EXPERIMENTAL DETAILS 

Calibration of R 

We determine the distance R between the two trap 
centers in the following way. A second aspheric lens fac- 
ing the one used to generate the traps coUimates back the 
trap beams; a / = 500 mm lens is then used to conjugate 
the trapping plane onto a CCD camera. We carefully 
calibrate the magnification of the system by measuring 
the displacement of the trap image when one changes the 
angle of incidence 9 of one of the trapping beams by a 
known amount. The resulting uncertainty in the abso- 
lute value of R is estimated to be below 5 %, due to the 
combined uncertainties in the determination of 9 and in 
the effective focal length {f^s = 10.0 ±0.1 mm) of the 
aspheric lens. Note that for our range of parameters, the 
dispersion in the interatomic distances due to the finite 
temperature of the atoms in the traps is negligible. For 
T — 50 /iK, the rms transverse size in our traps, which 
have a radial frequency uj^ 2n x 100 kHz, is about 
100 nm. The longitudinal size only contributes to the in- 
teratomic distance in second-order, and is thus negligible 
for the range of inter-trap distances considered here. 

Cancelation of stray electric fields 

Our setup was designed to minimize residual electric 
fields acting on the Rydberg atoms: (i) the metallic vac- 
uum chamber is relatively large, such that the glass sur- 
faces of viewports are about 10 cm away from the atoms, 
which are further shielded by a large metallic piece hold- 
ing the aspheric lenses; (ii) the surface of the aspheric 
lenses facing the atoms are coated with a ~ 200 nm-thick 
Indium-Tin Oxide (ITO) layer. Moreover, we have in- 
cluded a set of eight independent electrodes allowing us to 
apply electric fields pointing in any direction. Using the 
Stark spectroscopy of the |62D3/2,mj- = 3/2) state, we 
could determine that with all eight electrodes grounded, 
a residual field of ^ 150 mV/cm (pointing essentially 
along the x direction, see Fig. 1) was present. After ap- 
plying appropriate voltages to the eight electrodes, we 
obtain residual fields that we estimate, using Stark spec- 
troscopy of the more sensitive |82Z?3/2) state, to be well 
below ~ 5 mV/cm. 

THEORETICAL MODEL FOR THE TWO- ATOM 
DYNAMICS 

We model the dynamics of the two-atom system in the 
following way. Each atom i (i = 1, 2) is considered as a 



two-level atom with two internal states: the ground state 
\g)i and the Rydberg state |r)i, coupled by a near- 

resonant laser with Rabi frequency Vli. In the two- atom 
basis {I5.9), \rg), |rr)}, and using the rotating wave 
approximation, the hamiltonian of the system thus reads 

/ hvi2l2 nni/2 o \ 
hn2/2 -5 hn2/2 

hVLi/2 -5 Mli/2 • ^ > 
\ hVt2/2 hni/2 U^d^ -26 ) 

We omit the diagonal terms corresponding to interac- 
tion between two ground-state atoms (for a distance 
R—1 /xm the van der Waals interaction is in the 10~* Hz 
range, the magnetic dipole-dipole interaction in the mHz 
range) or between a ground-state and a Rydberg atom 
(at most in the 1 Hz range at i? = 1 /-tm), as they are 
completely negligible with respect to the other relevant 
energy scales of the problem. In practice, for all our data 
points, we have f2i — Q.2 to within 5%. 

Experimentally, we observe (especially for large Rabi 
frequencies f2.i) a significant damping of the Rabi oscil- 
lations, essentially due to ofi^-resonant spontaneous emis- 
sion via the bPi/2 intermediate state. To obtain the time 
dependence of the various populations of the two-atom 
system, we thus solve the equations of motion for the 
two-atom density matrix p 

p^hH,p]+C. (2) 
in 

We write the dissipative term as £ = £1 0ld2 +Idi ^£2, 
where 

is the dissipator for the two-level atom i (neglecting any 
dephasing, i.e. assuming T2 — 2Ti), expressed in the 
basis {\g)i, \r)i}, and Idj is the 2x2 identity matrix 
acting on atom j. We thus neglect, as far as dissipation 
is concerned, any cooperative effect between the atoms, 
such as super- or sub-radiance (this is legitimate as the 
atom separation R is much larger than the wavelength 
A ~ 795 nm of the 5S-5P transition which dominates the 
dissipation via spontaneous emission) [24]. 

NUMERICAL CALCULATION OF U^dw{R) 

To compute the interaction energy between two Ry- 
dberg atoms we diagonalized the interaction Hamilto- 
nian, considering only the leading, electric dipole interac- 
tion term. The diagonal terms of the Hamiltonian were 
computed by quantum defect theory, using the measure- 
ments reported in [SI]. The radial wavefunctions of the 
Rydberg states were numerically calculated using the 
Numerov algorithm [S2] and used to compute the off- 
diagonal terms of the interaction Hamiltonian. 



In order to deal with the very large number of Ry- 
dberg states involved, we restrict the Hamiltonian to a 
two-atom basis |ni/ij'imji, n2^2j2'Tij2) comprising only 
states which are close in energy (up to AE ^ ft, x 5 GHz) 
from the \nD^/2j''^j = 3/2, 711)3/2, — 3/2^ state, and 
with a principal quantum number difference An — 
I '^ — '^1,2! < 4. This corresponds to matrices with di- 
mensions up to ~ 10'^ X 10^. Doubling the size of the 
basis (by increasing AE and An) does not significantly 
(< 0.1%) alter the final results. 

The observed behavior is well described by the van der 
Waals interaction Cq/R^, and is not significantly affected 
by any level crossings, which are observed to occur at 
interatomic distances smaller than those experimentally 
explored in this paper. We obtain the values of the Cg 
coefficients by fitting the numerically obtained interac- 
tion energies U{R) at large distances 15 < ii < 20 /im; 
see Table below. We have also checked that the results 
are in good agreement with second-order perturbation 
theory [28]. 

In order to get a rough estimate of the precision in our 
determination of the Cg coefhcients, we have checked how 
numerical errors in the calculations of the radial matrix 
elements could impact their values. The uncertainties 
which are given in the table for jCgl^jj (see next section) 
were obtained by adding a random, uniformly-distributed 
relative error of ±0.5% to all radial matrix elements ap- 
pearing in the hamiltonian. We notice that the relative 
error is larger (^ f 0%) for n — 53, which is explained by 
the opposite signs of the leading perturbative terms in 
the sum giving Cg. 

DETERMINATION OF THE C'a COEFFICIENT 

Fitting the experimental data by |f7vdwl — ICel/R^ 
with |Cg| as an adjustable parameter, we obtain the val- 
ues in the table below (where the stated uncertainty is 
the Icr confidence interval on the estimate of |Cg| value). 

As can be seen, the experimental values agree well with 
the theoretical ones, especially when one considers the 
systematic uncertainty arising from the 5% uncertainty 
in the calibration of R. 



n 


|Cg|fit (GHz-/im6) 


|Cg|th. (GHz-^m^) 


53 


13.7 ± 1.2 


16.9 ± 1.7 


62 


730 ± 20 


766 ± 15 


82 


8500 ± 300 


8870 ± 150 
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